Highly Efficient CRISPR/Cas9 Mediated Gene Editing in Ocimum basilicum ‘FT Italiko’ to Induce Resistance to Peronospora belbahrii

Ocimum basilicum (sweet basil) is an economically important aromatic herb; in Italy, approximately 1000 ha of “Genovese-type” basil are grown annually in greenhouses and open fields and are subjected to Downy Mildew (DM) disease, caused by Peronospora belbahrii, leading to huge crop losses. Mutation of the Susceptibility (S) gene DMR6 (Downy Mildew Resistant 6) has been proven to confer a broad-spectrum resistance to DM. In this work, an effective Genome Editing (GE) approach mediated by CRISPR/Cas9 in O. basilicum ‘Italiko’, the élite cultivar used to produce “Pesto Genovese D.O.P”, was developed. A highly efficient genetic transformation method mediated by A. tumefaciens has been optimized from cotyledonary nodes, obtaining 82.2% of regenerated shoots, 84.6% of which resulted in Cas9+ plants. Eleven T0 lines presented different type of mutations in ObDMR6; 60% were indel frameshift mutations with knock-out of ObDMR6 of ‘FT Italiko’. Analysis of six T1 transgene-free seedlings revealed that the mutations of T0 plants were inherited and segregated. Based on infection trials conducted on T0 plants, clone 22B showed a very low percentage of disease incidence after 14 days post infection. The aromatic profile of all in vitro edited plants was also reported; all of them showed oxygenated monoterpenes as the major fraction.


Introduction
Ocimum basilicum, known as sweet basil, is one of the major evergreen multipurpose aromatic herbs. It belongs to the Lamiaceae family and is considered allotetraploid (2n = 4x = 48) [1,2]. It is an economically important crop due to its chemical properties used in pharmaceutical, cosmetic, and culinary fields [3,4]. Leaves were extensively studied for its secondary metabolite content, which includes polyphenols such as flavonoids and anthocyanins [5], involved in important biological activities as antioxidant, anticancer, antimicrobial, anti-inflammatory, and insecticide. Furthermore, basil is appreciated for culinary flavoring and aromatherapy [6].

Identification of ObDMR6 in 'FT Italiko' and Protein Alignment
The sequence ObDMR6 in 'FT Italiko', obtained by BLASTN analysis, reached a max score alignment (371) and with E value of 4 × e −99 and showed sequence homology with AtDMR6. The translation of transcript was confirmed with BLASTX search in nonredundant protein database, showing 88.10% of identity with 2-oxoglutarate and oxygenase superfamily protein of Perilla frutescens var. frutescens (GenBank accession: KAH6789542.1), 87.50% of identity with downy mildew resistance-like protein 6 v3 and v1 of O. basilicum cv 'Genoveser' (GenBank accession: QWT44769.1 and QWT44767.1, respectively).
Whole genomic sequence of 1260 bp of ObDMR6 in 'FT Italiko' was isolated and sequenced. The alignment with AtDMR6 coding sequence suggested the presence and the position of four exons and three introns ( Figure 1). channel assuming the possible tetraploid content (2n = 48) (S1 file). This was also confirmed by the mean of fluorescence ratio between sample and reference calculated (approx. 1.9).

Identification of ObDMR6 in 'FT Italiko' and Protein Alignment
The sequence ObDMR6 in 'FT Italiko', obtained by BLASTN analysis, reached a max score alignment (371) and with E value of 4e-99 and showed sequence homology with AtDMR6. The translation of transcript was confirmed with BLASTX search in nonredundant protein database, showing 88.10% of identity with 2-oxoglutarate and oxygenase superfamily protein of Perilla frutescens var. frutescens (GenBank accession: KAH6789542.1), 87.50% of identity with downy mildew resistance-like protein 6 v3 and v1 of O. basilicum cv 'Genoveser' (GenBank accession: QWT44769.1 and QWT44767.1, respectively).
Whole genomic sequence of 1260 bp of ObDMR6 in 'FT Italiko' was isolated and sequenced. The alignment with AtDMR6 coding sequence suggested the presence and the position of four exons and three introns ( Figure 1). The resulting ObDMR6 cds was predicted to have an ORF of 1008 bp encoding 336 amino acids. The ObDMR6 'Ft Italiko' sequence was deposited on NCBI GenBank with the accession number MT319764 and released on 30 September 2021.
The CoGeBlast analysis has located the blast hits in six chromosome scaffolds (262, 7503, 4335, 8208, 235, 9553), discriminating almost six versions of ObDMR6, as described for sweet basil cv 'Genoveser' by [38]. These sequences show different SNPs against ObDMR6 'FT Italiko' and among them, only scaffold 7503 and scaffold 262 presented high percentage of identities with query sequence, respectively, of 95% and 93%.

Editing in Hairy Roots
To evaluate the efficiency of the CRISPR/Cas9 system in basil, genetic transformation experiments mediated by Agrobacterium rhizogenes were performed. A high induction of Hairy Roots (HR) was obtained from all explants tested (cotyledons, leaves, in vitro stems) (S5 file). Twenty-nine out of thirty HR were found to be positive in the PCR analysis for the presence of Cas9, confirming a high rate of co-transformation equal to 96% (introduction of the T-DNA genes of the Ri plasmid and the pDirect 22C plasmid).
Sequencing of the 396 bp fragment comprising sgRNAs target site ( Figure 2) of two Cas9 + HR lines, allowed us to identify frameshift mutations caused by the deletion of a nucleotide (C), in position 5 nt upstream of the PAM (Figure 3). The resulting ObDMR6 cds was predicted to have an ORF of 1008 bp encoding 336 amino acids. The ObDMR6 'Ft Italiko' sequence was deposited on NCBI GenBank with the accession number MT319764 and released on 30 September 2021.
The CoGeBlast analysis has located the blast hits in six chromosome scaffolds (262, 7503, 4335, 8208, 235, 9553), discriminating almost six versions of ObDMR6, as described for sweet basil cv 'Genoveser' by [38]. These sequences show different SNPs against ObDMR6 'FT Italiko' and among them, only scaffold 7503 and scaffold 262 presented high percentage of identities with query sequence, respectively, of 95% and 93%.

Editing in Hairy Roots
To evaluate the efficiency of the CRISPR/Cas9 system in basil, genetic transformation experiments mediated by Agrobacterium rhizogenes were performed. A high induction of Hairy Roots (HR) was obtained from all explants tested (cotyledons, leaves, in vitro stems) (S5 file). Twenty-nine out of thirty HR were found to be positive in the PCR analysis for the presence of Cas9, confirming a high rate of co-transformation equal to 96% (introduction of the T-DNA genes of the Ri plasmid and the pDirect 22C plasmid).
Sequencing of the 396 bp fragment comprising sgRNAs target site ( Figure 2) of two Cas9 + HR lines, allowed us to identify frameshift mutations caused by the deletion of a nucleotide (C), in position 5 nt upstream of the PAM (Figure 3).

Production of T0 Edited Plants by Agrobacterium Tumefaciens-Mediated Transformation
After four weeks of culture, it was possible to detect shoots that directly emerged from 150 CNs, without any visible callus formation ( Figure 4A,B). Eighty-two percent of regeneration efficiency, calculated as the ratio between number of regenerated CNs on total number of CNs, was recorded on kanamycin containing medium, with a mean of 2.6 of regenerated shoots for each explant. CNs, excided from 4-week-old plantlets generated transgenic plants expressing sgRNA(s) and Cas9.

Production of T0 Edited Plants by Agrobacterium Tumefaciens-Mediated Transformation
After four weeks of culture, it was possible to detect shoots that directly emerged from 150 CNs, without any visible callus formation ( Figure 4A,B). Eighty-two percent of regeneration efficiency, calculated as the ratio between number of regenerated CNs on total number of CNs, was recorded on kanamycin containing medium, with a mean of 2.6 of regenerated shoots for each explant. CNs, excided from 4-week-old plantlets generated transgenic plants expressing sgRNA(s) and Cas9. Different clones were obtained by shoot multiplication ( Figure 4C). After 30 days, all singularized plants were able to develop an efficient root system ( Figure 4D). Interestingly, some in vitro clones bloomed ( Figure 4E). All clones transferred in insectproof greenhouse were acclimatized and three T0 lines (21D, 22B and 45B) flowered and produced seeds ( Figure 4F). Twenty-six T0 plants were tested for Cas9 presence by PCR, and 84.6% resulted to be positive for the integration of the transgene.  Different clones were obtained by shoot multiplication ( Figure 4C). After 30 days, all singularized plants were able to develop an efficient root system ( Figure 4D). Interestingly, some in vitro clones bloomed ( Figure 4E). All clones transferred in insect-proof greenhouse were acclimatized and three T0 lines (21D, 22B and 45B) flowered and produced seeds ( Figure 4F). Twenty-six T0 plants were tested for Cas9 presence by PCR, and 84.6% resulted to be positive for the integration of the transgene.
To identify the mutations in 21 ObDMR6 transgenic plants, a 396 bp fragment comprising both sgRNA target site ( Figure 2) was amplified by PCR and subsequently Sanger sequenced for each samples analyzed.
The sequence alignment of 396 bp fragments among fourteen WT individual clones of ObDMR6 showed four different variants (S6 file). The SNP (A/T) at 11 bases upstream the PAM leads to amino acid change in the sequence: tyrosine (Y)/phenylalanine (F).
Different type of mutations in ObDMR6 were identified in 86% of T0 lines transformed with pDirect22c_sg442-462 ( Figure 4).
Seven different mutation patterns were observed in the T0 plants. Seven out of fifteen T0 mutants had nucleotide deletions around the target site, while two mutants show insertions of one and two bases within target site.
More than sixty percent (9 out of 15) of the T0 mutants were confirmed to have frameshift mutations that cause the knock-out of ObDMR6 gene and the premature termination of gene product. Five out of seventeen of the mutant plants show one or more base substitutions leading to amino acid sequence changes (S4 File).

Cas9-Free ObDMR6 Mutants in T1 Generation
Six T1 transgene-free seedlings (derived from 21D and 22B lines), screened by PCR, having eliminated the Cas9 gene by segregation, were analyzed for the ObDMR6 mutation derived from editing and compared to WT, by direct Sanger sequencing of the 396 bp fragment (almost three fragments for each plant were analyzed). The chromatograms generated for each plant were compared to WT and among T0 and T1 of the same line. Plants showing chromatograms with overlapping spikes (two nucleotides simultaneously) were considered heterozygous.
The chromatograms obtained from each T1 plants reveal that the mutations of T0 plants were inherited and segregated. Plants of line 22B had the deletion of a nucleotide (C), in position 5 nt upstream of the PAM, leading to the creation of a premature termination (stop) codon. For the T1 plants analyzed, the mutation appears to be homozygous, although the chromatogram of plant T1_22B_6 also shows a heterozygous single nucleotide polymorphism (SNP) that is anyway outside of the editing region: two peaks of equal heights but different colors (A/T) (as highlighted with ** in Figure 5). Plants of line 21D had the substitution of two nucleotides CG in place of TT, in position 7 nt upstream from the PAM, resulting in amino acid change in the sequence: alanine (A) instead of valine (V). Plant T1_21D_1 inherited the same mutation of T0_21D and displayed further nucleotide changes upstream of the PAM (CG in place of TG). Furthermore, it showed the same heterozygous SNP of plant T1_22B_6 (also denoted by **).
Seven different mutation patterns were observed in the T0 plants. Seven out of fifteen T0 mutants had nucleotide deletions around the target site, while two mutants show insertions of one and two bases within target site.
More than sixty percent (9 out of 15) of the T0 mutants were confirmed to have frameshift mutations that cause the knock-out of ObDMR6 gene and the premature termination of gene product. Five out of seventeen of the mutant plants show one or more base substitutions leading to amino acid sequence changes (S4 File).

Cas9-Free ObDMR6 Mutants in T1 Generation
Six T1 transgene-free seedlings (derived from 21D and 22B lines), screened by PCR, having eliminated the Cas9 gene by segregation, were analyzed for the ObDMR6 mutation derived from editing and compared to WT, by direct Sanger sequencing of the 396 bp fragment (almost three fragments for each plant were analyzed). The chromatograms generated for each plant were compared to WT and among T0 and T1 of the same line. Plants showing chromatograms with overlapping spikes (two nucleotides simultaneously) were considered heterozygous.
The chromatograms obtained from each T1 plants reveal that the mutations of T0 plants were inherited and segregated. Plants of line 22B had the deletion of a nucleotide (C), in position 5 nt upstream of the PAM, leading to the creation of a premature termination (stop) codon. For the T1 plants analyzed, the mutation appears to be homozygous, although the chromatogram of plant T1_22B_6 also shows a heterozygous single nucleotide polymorphism (SNP) that is anyway outside of the editing region: two peaks of equal heights but different colors (A/T) (as highlighted with ** in Figure 5). Plants of line 21D had the substitution of two nucleotides CG in place of TT, in position 7 nt upstream from the PAM, resulting in amino acid change in the sequence: alanine (A) instead of valine (V). Plant T1_21D_1 inherited the same mutation of T0_21D and displayed further nucleotide changes upstream of the PAM (CG in place of TG). Furthermore, it showed the same heterozygous SNP of plant T1_22B_6 (also denoted by **).

Morphological Phenotyping of T0 and T1 Plants
The morphological traits observed in T0 plants are reported in Table 1 and S7 File. For each clone and WT were evaluated almost three plants. The clones 21D, 47C, and 45B showed strongly decrease in plant height, internode number and leaf length. The shorter internode length and curling up margin of some leaves characterized 21D and 47C clones. The clones 17A2 and 32A differed from WT ones with respect to their height. Their leaves showed a normal size, but few assuming a shape with curling up margin. Plants of 22B clone are similar to WT plants, except for the morphology of few leaves. Flowering was normal in 22B, 21D, and 45B, but absent in the other three clones.

P. belbahrii Infection Trial
As reported in Tables 3 and 4, the values of Disease Incidence and Disease Severity increased in percentage from 7 to 14 dpi for all the clones, reaching the maximum value (100%) for WT, 17A, 47C, and 21D. Even none of the clones showed infection symptoms at the end of the trial, clones 47B and 22B displayed a medium level of tolerance to the infection. Particularly, at 7 dpi, clone 22B did not show any infection symptom (0% of DI) and very low percentage of incidence and severity were detected at 10 and 14 dpi, in respect to the others clones and WT, confirming that this clone could be very promising. Visual symptoms on leaves of the clones 47B and 22C, showing the absence of spots and or chlorosis on the underside of the leaves and compared to the WT are reported in Figure 6. Table 3. Percentage of infected leaves (Disease Incidence) at 7-, 10-and 14-days post inoculation (dpi) (different letters indicate significant difference among clones within dpi (p ≤ 0.05, Tukey HSD test).

Volatile Profile Analysis
The complete chemical compositions of the headspaces (HSs) of the analyses O. basilicum clones are reported in Table 5. A total of 36 compounds were identified, accounting for 99.6-100.0% of the whole compositions. The class of terpenes was the most represented, followed by phenylpropanoids. Among the former class, monoterpenes were detected in good amounts in both their hydrocarbons and oxygenated forms, even though the latter was sharply predominant, accounting for 33.8-64.5%. The control sample was

Volatile Profile Analysis
The complete chemical compositions of the headspaces (HSs) of the analyses O. basilicum clones are reported in Table 5. A total of 36 compounds were identified, accounting for 99.6-100.0% of the whole compositions. The class of terpenes was the most represented, followed by phenylpropanoids. Among the former class, monoterpenes were detected in good amounts in both their hydrocarbons and oxygenated forms, even though the latter was sharply predominant, accounting for 33.8-64.5%. The control sample was characterized by the highest amount, followed by the pairs 21D (47.4%) and 47C (47.1%) and 17A2 (39.4%) and 45B (36.3%), each with similar contents of oxygenated monoterpenes. Conversely, the sample 22B (33.8%) showed the lowest relative abundance of these secondary metabolites. 1,8-Cineole was the most important compound detected in all the samples, whose content ranged from 29.9 to 44.1%. Linalool, instead, was found in percentages above 15% in the control HS, while all the other samples showed narrowed contents, always below 5%. Monoterpene hydrocarbons, mainly represented by (E)-β-ocimene (3.2-11.2%) e βpinene (2.2-4.9%), were also found in not negligible amounts, as they ranged from 16.6 to 24.6% of the HSs. The highest content of this class was detected in the 21D HS, while the significantly lowest in those of 45B and CTR.
Among terpenes, also sesquiterpene hydrocarbons represented an important chemical class, whose content was comprised between 8.1 and 21.4%. These compounds were minimum in the control sample, while their content was greatest in 45B. All the other samples, instead, did not show significant differences. Among sesquiterpene hydrocarbons, trans-α-bergamotene was certainly the most relevant one.
Finally, the class of phenylpropanoids was detected in percentages from 10.7 in CTR up to 29.4% in 22B. This class was represented only by eugenol and its methylated derivative. The former compound was higher in all the analyzed clones, with the only exception of 45B, in which, instead methyl eugenol prevailed.
The complete chemical composition of the HSs was subjected to multivariate statistical analysis with the Hierarchical Cluster Analysis (HCA) and the Principal Component Analysis (PCA) methods whose dendrogram and plots are reported in Figures 7 and 8, respectively.
Among terpenes, also sesquiterpene hydrocarbons represented an important chemical class, whose content was comprised between 8.1 and 21.4%. These compounds were minimum in the control sample, while their content was greatest in 45B. All the other samples, instead, did not show significant differences. Among sesquiterpene hydrocarbons, trans-α-bergamotene was certainly the most relevant one.
Finally, the class of phenylpropanoids was detected in percentages from 10.7 in CTR up to 29.4% in 22B. This class was represented only by eugenol and its methylated derivative. The former compound was higher in all the analyzed clones, with the only exception of 45B, in which, instead methyl eugenol prevailed.
The complete chemical composition of the HSs was subjected to multivariate statistical analysis with the Hierarchical Cluster Analysis (HCA) and the Principal Component Analysis (PCA) methods whose dendrogram and plots are reported in Figures 7 and 8, respectively. The dendrogram of the HCA was characterized by two main clusters, the red and the green ones. The former comprised the samples CTR, 21D, and 47C, of which the last two showed a great similarity, gathering in the same subgroup. The same was found for 17A2 and 22B, which along with 45B, formed the green cluster. A similar partition was also evidenced by the score plot of the PCA (Figure 8a), in which the samples belonging to the green group of the dendrogram were plotted in the left quadrants (PC1 < 0), while those of the red cluster in the right ones (PC1 > 0). Regarding the former group, the positioning of the samples 22B and 17A2 was, perhaps, determined by the greatest content of eugenol, while that of 45B by methyl eugenol. Conversely, the samples of the red group were probably plotted in the right quadrants on behalf of the 1,8-cineole vector. However, the positioning of the control sample in the rightmost area of the upper quadrant (PC2 > 0) was attributable to the linalool vector, whose content was significantly greatest in its HS. Both statistical analyses evidenced a higher similarity of the control sample to the clones 47C and 21D, in terms of spontaneous volatile emission, as they are all characterized by both a greater content of oxygenated monoterpenes and a lower amount of phenylpropanoids.

Discussion
Climate change is causing rapid changes in the intensity and spread of biotic stresses, such as old and new parasites [40]. It is, therefore, essential to have adequate technologies to develop new plant genotypes endowed with durable and broad-spectrum resistance. Traditional breeding does not always succeed in producing new genotypes capable of coping with rapid changes. Gene-editing technology makes possible to modify target DNA in a specific site, in a precise and fast way, inducing specific genetic variations, but maintaining the genetic identity [41] of elite cultivars. The inactivation of susceptibility genes required host-pathogen recognition or supporting infection process at different step and has allowed us to obtain plants resistant to different pathogens, in several economically important crops. Furthermore, the CRISPR/Cas9 technology makes it possible to generate transgene-free gene-edited plants resistance, which are better accepted by public opinion and do not have to fall back into European GMO legislation [42].
In this work, an effective Genome Editing approach mediated by CRISPR/Cas9 in O. The dendrogram of the HCA was characterized by two main clusters, the red and the green ones. The former comprised the samples CTR, 21D, and 47C, of which the last two showed a great similarity, gathering in the same subgroup. The same was found for 17A2 and 22B, which along with 45B, formed the green cluster.
A similar partition was also evidenced by the score plot of the PCA (Figure 8a), in which the samples belonging to the green group of the dendrogram were plotted in the left quadrants (PC1 < 0), while those of the red cluster in the right ones (PC1 > 0). Regarding the former group, the positioning of the samples 22B and 17A2 was, perhaps, determined by the greatest content of eugenol, while that of 45B by methyl eugenol. Conversely, the samples of the red group were probably plotted in the right quadrants on behalf of the 1,8-cineole vector. However, the positioning of the control sample in the rightmost area of the upper quadrant (PC2 > 0) was attributable to the linalool vector, whose content was significantly greatest in its HS.
Both statistical analyses evidenced a higher similarity of the control sample to the clones 47C and 21D, in terms of spontaneous volatile emission, as they are all char-acterized by both a greater content of oxygenated monoterpenes and a lower amount of phenylpropanoids.

Discussion
Climate change is causing rapid changes in the intensity and spread of biotic stresses, such as old and new parasites [40]. It is, therefore, essential to have adequate technologies to develop new plant genotypes endowed with durable and broad-spectrum resistance. Traditional breeding does not always succeed in producing new genotypes capable of coping with rapid changes. Gene-editing technology makes possible to modify target DNA in a specific site, in a precise and fast way, inducing specific genetic variations, but maintaining the genetic identity [41] of elite cultivars. The inactivation of susceptibility genes required host-pathogen recognition or supporting infection process at different step and has allowed us to obtain plants resistant to different pathogens, in several economically important crops. Furthermore, the CRISPR/Cas9 technology makes it possible to generate transgene-free gene-edited plants resistance, which are better accepted by public opinion and do not have to fall back into European GMO legislation [42].
In this work, an effective Genome Editing approach mediated by CRISPR/Cas9 in O. basilicum 'Italiko', the elite cultivar used to produce "Pesto Genovese D.O.P", was developed.
The regeneration of whole plant from modified cells and tissues, through in vitro culture techniques, represents the main bottleneck for the application of Genome Editing in different plant species. A high rate of regenerated shoots is a necessary requirement for efficient genetic transformation. High percentage of sprouts regenerated by the CNs in O. basilicum cv 'Italiko' was recorded, confirming that cotyledonary node-based direct regeneration method is effective and suitable for genetic transformation mediated by A. tumefaciens [43]. Furthermore, direct shoot organogenesis, without callus proliferation, shortens the time to obtain regenerated plants and decreases the risk of somaclonal variation [44], facilitating the selection of transgenic plants.
A successful transformation protocol mediated by A. rhizogenes, carrying pDirect 22C binary vector and Ri plasmid, allowed us to obtain HRs with a high rate of cotransformation. Several studies reported the effective application of A. rhizogenes system in genome editing [45,46]. HRs, thanks to their prolific and quick growth rate and being easy to maintain for long time, represent a robust tool to test the efficiency of CRISPR/Cas9 system. Furthermore, the recalcitrant-to-regeneration plants could display after genome editing a good shoot regeneration potential [47].
Several studies conducted in different species of plants and crops (tomato, potato, barley, banana, citrus, apple, grapevine, basil) have shown that the mutation of S gene DMR6, discovered and characterized in A. thaliana [24], conferred resistance to oomycete and other pathogens. Hasley and colleagues [38], targeting ObDMR6 in sweet basil 'Genoveser' by CRISPR/Cas9, found enhanced DM resistance with reduced production of sporangia and pathogen accumulation.
The sequence of ObDMR6 from 'FT Italiko' was deposited on NCBI GenBank at the same time with the release of sequences of sweet basil DMR6 cv Genoveser, (Enza Zaden) [38]. The comparison with the latter and with the O. basilicum genome sequences recently released [48] indicated the presence of six variable copies of the DMR6 gene in basil, being the 'FT Italiko' sequence used in the present study closer to copies present in genomic scaffolds 7503 and 262. O. basilicum presents high interspecific variability even among cvs, in terms of morphological and chemical traits (aromatic profile) and susceptibility or resistance gene against pathogens [1,49].
The two overlapping sgRNA target sequences were designed on the second exon of the gene and were found to be conserved in all six chromosome scaffolds identified by CoGeBlast tool. It is reported that the overlapping sgRNAs (that share some base pairs) create small nucleotide deletions in target site with higher efficiency [50]. The presence of a SNP (A/T) 11 bases upstream one of the two PAMs in a sample of 14 wild type 'FT Italiko' plants should not reduce Cas9 ability to bind and cleave its target [51]. Indeed, different type of mutations in ObDMR6 were identified in 86% of tested T0 transgenic plants; the 60% are indel frameshift mutations with knock-out of ObDMR6 'FT Italiko' and essentially two types of deletions (−1 and −2) and insertions (+1 and +2). As reported in two recent manuscripts [38,52], although basil is tetraploid, it has been possible obtain homozygous mutant plants free of the transgene at T1 generation also in O. basilicum cv 'Italiko'. Chromatograms generated for each plant of two T1 lines (21D and 22B) clearly showed that mutations were inherited from their parents and were homozygous. Differences between the sequences of T0 plants and T1 progenies may be due to Cas9 activity between the time of tissue sampling of T0 plants and gene segregation following meiosis. Additionally, ref. [53] reported different mutation pattern between T1 plants and their own parent plants.
The evaluation of some morphological traits measured on acclimatized in primary transgenic plants (T0), showed differences among clones and WT, probably as a consequences of in vitro culture and/or DMR6 gene editing. Only the 22B clone appears similar to WT, for plant and leaf traits. Few papers described the phenotype of DMR6 mutants in different crops. In potato, an affection in growth was reported in phenotypes for StDMR6-2 mutants, with significantly lower plant height and shorter internodes [34]; while in tomato [31], SlDMR6-1 mutants displayed indistinguishable phenotype from WT plants under laboratory conditions. No phenotype description in sweet basil has been reported for ObDMR6 mutants [38].
To the best of our knowledge, this is the first study reporting the chemical composition of the spontaneous volatile emission of in vitro plants of O. basilicum 'FT Italiko' and the first concerning the aroma of edited sweet basil plants. The aroma composition represented one of the most important aspects of O. basilicum cultivars used for the preparation of pesto sauce, as it is responsible for its organoleptic properties.
The main chemical compounds detected in the control sample partially corroborate the findings of [54] who reported a prevalence of linalool, instead of 1,8-cineole in the volatile composition of O. basilicum cv. Genovese. In [55], oxygenated monoterpenes are the major secondary metabolites found in vitro plants, with abundance of linalool and 1,8 cineole depending on the different culture media used. According to our findings, gene editing was able to significantly influence the aroma profile of basil plants, creating two groups of samples: the former including the clones 21D and 47C, more similar to the control sample, but with a marked lower content of linalool on behalf of eugenol, and the second, including the clones 17A2, 22B, and 45B, in which the synthesis of methyl eugenol was strongly enhanced. Methyl eugenol is a common phenylpropanoid found in several plant species, whose release is seemingly related to the plant defense against pathogens and herbivore insects [56]; this information could be related to the fact that plants 45B and 22B are less susceptible to the Peronospora infection.
In the experimental trial with P. belbahrii fresh spore solution, under permissive conditions, T0 22B and 21D plants have shown an enhanced DM resistance than WT, confirming the involvement of the DMR6 gene in the interaction and recognition with its pathogen and that its inactivation improves resistance to DM. Quite significant differences are observed among Cas9+ plants in disease resistance, morphological traits and aromatic profile are probably due to intra-variety variability of commercial seeds as starting material for obtaining CNs and to independent origin of each regenerant plant from CNs following genetic transformation. Furthermore, other factors should be considered, as non-predictable integration of non-T-DNA portions of the vector into the plant genome [57] due incorrect recognition of borders by Agrobacterium and off-target effects on the genome that were predicted only in silico. In addition, multiple members of the DMR6-like gene family, recently identified in sweet basil [38,58], may contribute to showing varying degrees of silencing and partially restore the susceptibility to DM. In the future, the progenies derived from these promising plants will have to be tested in field trials to evaluate their potential in comparison to current and original cv and subsequently used in breeding programs to create new lines of elite cv sweet basil resistant to DM.

Flow Cytometry
Ploidy level of O. basilicum 'FT Italiko' was estimated using cytometry. Cell nuclei were isolated from upper young leaves of basil plant, together with reference plant (O. tenuiflorum = O. sanctum, 2C = 0.72 pg), by chopping leaves with a sharp razor blade in the Petri dish in 500 µL of nuclei extraction buffer (CyStain-UV precise P, Sysmex-Partec). To reduce the negative effect of cytosolic and phenolic compounds, 1% of PVP-40 (polyvinylpyrrolidone-40; Sigma PVP40) was added to the buffer [59]. Then, the crude suspension was filtered through a 30 µm CellTrics disposable nylon filter and was kept on ice for 10 min. For samples coloring, 1.6 mL of DAPI staining solution (4 ,6-diamino-2-phenylindole) was added, according to the manufacturer's instructions. The fluorescence intensity was measured by CyFlow Ploidy Analyser (Sysmex-Partec), equipped with a UV LED chip (365 nm) for DAPI excitation. The flow rate varied between 20 and 50 events per second, and about 6000 nuclei per sample were measured. Histograms were acquired using CyView software (Sysmex-Partec). Reference peak was gated manually, and coefficient of variation [60] and maximum peak height were recorded. The mean of fluorescence ratio between sample and reference was used to estimate the DNA ploidy levels. Analyses were performed twice for each sample (5 leaf samples of O. basilicum 'FT Italiko' and 5 leaves of O. tenuiflorum plants).

Regeneration Protocol Via Direct Organogenesis and Agrobacterium Tumefaciens Mediated Transformation
Seeds of Ocimum basilicum L. 'FT Italiko' (La Semiorto Sementi ® ) were sterilized and germinated according to [61] and cotyledonary nodes (CNs) were used as starting material for plant transformation and direct shoot formation as suggested by [62] on regeneration medium (4.3 g/L MS "shoot multiplication medium", 3% sucrose, 3.7 mg/L TDZ, 8 g/L agar, pH 5.8).
A. tumefaciens was cultured in 5 mL of liquid TYNG medium (10 g/L tryptone, 5 g/L yeast extract, 5 g/L NaCl, and 0.5 g/L MgSO 4 ·7H 2 O, pH 7.5); supplemented with rifampicin 50 mg/L, kanamycin 100 mg/L, and carbenicillin 100 mg/L; and incubated at 28 • C for 24 h at 120 r.p.m. The culture was incubated overnight under the same conditions and then centrifuged at 2000× g for 15 min at room temperature. Pellet was resuspended in MS liquid medium up to OD600 of 0.25-0.3; acetosyringone was added to a final concentration of 100 µM before plant inoculation.
One-hundred fifty CNs were wounded, immersed for 30 min in suspension of A. tumefaciens cells, blotted dry and then placed on regeneration medium. After three days of co-cultivation, CNs were transferred to a regeneration medium supplemented by cefotaxime 100 mg/L and kanamycin 30 mg/L and sub-cultured every two weeks.
Regenerated shoots developed from the CNs were excised and cultivated on a multiplication medium (4.3 g/L MS, 3% sucrose, 0.2 mg/L 6-Benzyladenine (BA), 8 g/L agar, pH 5.7) with addition of cefotaxime 100 mg/L and kanamycin 50 mg/L. Each shoot was considered as a single clone. All plant materials were maintained in vitro at 23 ± 1 • C with a 16 h photoperiod of light (PPFD 30 µE m −2 s −1 ).

Agrobacterium Rhizogenes Mediated Transformation
A. rhizogenes strain LB9402 was cultured in 5 mL Nutrient Broth (NB) medium supplemented with kanamycin 100 mg/L according to the protocol described above.
Different explants (leaves, internodes, hypocotyls, petioles) were excised from seedlings, wounded, and immersed in bacterial suspension for 20 min at room temperature and then dried with sterile filter paper. The explants were placed on MS solid co-cultivation medium for 3 days at 24 • C in dark condition and subsequently placed on the same substrate with the addition of cefotaxime 100 mg/L and kanamycin 50 mg/L. After two weeks, 30 hairy roots (of 3 cm) were excised and kept separate as independent lines; they were sub-cultured every 15 days onto MS medium enriched with antibiotics.

Identification of ObDMR6 in 'FT Italiko'
In order to identify ObDMR6 in 'FT Italiko', AtDMR6 coding sequence was used as a query (GenBank accession: NM_122361) against the Transcriptome Shotgun Assembly Sequence Database of Ocimum tenuiflorum [63]. The sequence with the highest alignment score and with 4 × e −99 value was checked in NCBI (https://www.ncbi.nlm.nih.gov/, accessed on 1 March 2021) non-redundant protein database (nr) using a BLASTX search, to confirm the translation.
Genomic DNA was extracted from fresh 'FT Italiko' basil leaves, using Dneasy ® Plant Mini Kit (Qiagen, Hilden, Germany) according to manufacturer's protocol and used as template to perform PCR with Q5 High-Fidelity DNA Polymerase (NEB, Ipswich, MA, USA) following these conditions: initial denaturation at 98 • C for 2 min; followed by 35 cycles at 98 • C for 30 s, 55 • C for 30 s and 72 • C for 30 s, with a final extension at 72 • C for 4 min.
The PCR product (1260 bp) was cloned in pJET1.2/blunt vector (ThermoFisher ® , Scientific Inc., Waltham, MA, U.S.) and subjected to Sanger sequencing. The genomic nucleotide sequence was aligned with AtDMR6 cds to detect exons and introns position; the resulting ObDMR6 cds was translated using Expasy tool (www.expasy.org, accessed on 15 April 2021) and the amino acid sequence was searched in nr protein database using BLASTP to assess its homology. The amino acid sequence alignment was also generated using CLUSTALW tool of Bioedit 7.2 software (https://bioedit.software.informer.com, first access on 15 April 2021).

Protein Alignment
The alignment of six ObDMR6 amino acid variants (QWT44767.

Identification of CRISPR/Cas9 Targets and Vector Construction
Two sgRNA target sequences were identified in the coding region by CRISPR Direct design tool (https://crispr.dbcls.jp/, accessed on 15 April 2021) [65] and selected for editing 'FT Italiko' ObDMR6: sgRNA442 (target sequence CAGGCTTCATTGCTACCCCT) and sgRNA462 (target sequence TGGAGAATTATGTTCCTGAA) directed to the second exon of the gene). The sgRNA target sequences are specific for scaffold 7503 but presents a SNP in position 10 of sgRNA462 target for the scaffold 262 (CoGe database).
The vector was sequenced in order to verify the correct insertion of the cassette, then it was introduced into A. tumefaciens strain AGL1 [67] for O. basilicum 'FT Italiko' transformation, according to the protocol described by [68]. pDirect22c_sg442-462 was also introduced in A. rhizogenes strain LB9402 to evaluate the efficiency of genome editing in the élite cv of O. basilicum hairy roots system.

Screening of Transgenic Plants and HRs for Target ObDMR6 Mutations
Genomic DNA was extracted from 50 mg of HRs and from 100 mg leaf tissue of newly regenerated basil plantlets (T0) and wild-type plants, using DNeasy ® Plant mini-Kit, (QIAGEN, Hilden, Germany). In order to detect transgene integration, PCR was performed using primers Cas9_Fw and Cas9_Rv (5 -GTGCAGACCTACAACCAGCT and 5 -CTGAGAGTGGAGCCTTGGTG).
PCR was performed using Q5 high fidelity DNA polymerase (NEB) and PCR products were cloned directly in pJET1.2/blunt vector (ThermoFisher ® Scientific Inc., Waltham, Massachusetts, U.S.), followed by Sanger sequencing using primers provided by the cloning kit. At least three colonies for each clone have been sequenced.
Nucleotide sequence alignment was generated using CLUSTALW tool of Bioedit 7.2 software (https://bioedit.software.informer.com, access on 15 December 2021) and the detection of insertion-deletion mutations (indels) was performed.

Phenotypic Evaluation of T0 and T1 Mutant Plants
Newly regenerated transgenic T0 seedlings were in vitro multiplied and rooted. Phenotypic consequences of transformation were observed in primary transgenic plants (T0), which were potted and acclimatized in insect-proof greenhouse. Phenotypic traits, such as plant height, internode number and length, leaf morphology, length, and width, were evaluated and measured at the vegetative stage, on almost three plants of each clone.
Selfing bags were mounted on flower stalks of T0 plants, at the beginning of flowering time, to avoid cross-pollination and mature seeds were harvested.
Fifty seeds obtained from self-pollination of T0 plants 21D and 22B were potted and grown in greenhouse in comparison to wild type plants ('FT Italiko') at the same age. The aforementioned phenotypic traits were also recorded and statistically evaluated, using Student's t-test, in fifteen T1 plants, arising from the two clones.

Evaluation of In Vivo T0 Plants Infection
Preliminary infection trial with P. belbahrii on T0 plants to evaluate the acquired resistance through editing of the ObDMR6 was performed under controlled conditions at CeRSAA (Albenga, Italy).
Suspension of spores of known concentration (10 6 /mL) obtained from affected fresh leaves showing evident signs of sporulation was used as inoculum and it was sprayed directly on in vivo T0 plants. Symptoms were recorded at 7-, 10-and 14-days post inoculation (dpi). Disease Incidence (DI: % of infected leaves) and Disease Severity (DS: % of infected leaf surface) were calculated. Infection (sporulation) was evaluated in Selected  The headspace solid phase microextraction (HS-SPME) was used for the analyses of the spontaneous volatile emission of the O. basilicum plants. All the samples from in vitro T0 and WT plants (5 leaves each) were put in a 100 mL glass beaker covered with aluminum foil and left to equilibrate for at least 30 min. Then, the headspaces (HSs) were sampled with a Supelco PDMS fiber (100 µm) (Supelco analytical, Bellefonte, PA, USA), preconditioned according to the manufacturer's guidelines, for 10 min at room temperature, and once the exposure time was finished the fiber was withdrawn into the needle and suddenly injected into the GC-MS equipment. The analyses were accomplished in triplicates.

Gas Chromatography-Mass Spectrometry Analyses
The gas chromatography-electron impact mass spectrometry (GC-EIMS) analyses were performed with an Agilent 7890B gas chromatograph (Agilent Technologies Inc., Santa Clara, CA, USA) equipped with an Agilent HP-5MS capillary column (30 m × 0.25 mm; coating thickness 0.25 µm) and an Agilent 5977B single quadrupole mass detector using the following analytical conditions: oven temperature ramp rising from 60 to 240 • C at 3 • C/min; injector temperature at 250 • C; transfer line temperature at 240 • C; carrier gas helium, with a 1 mL/min flow. The acquisition parameters were set as follows: full scan; scan range: 30-300 m/z; scan time: 1.0 s. The identification of the constituents was based on a comparison of the retention times with those of the authentic samples, comparing their linear retention indices relative to the series C8-C27 of n-hydrocarbons. Moreover, computer matching was used against commercial (NIST 14 and ADAMS 2007) and laboratory-developed mass spectra libraries built up from pure substances and components of commercial essential oils of known composition and MS literature data [69].

Statistical Analysis
Multivariate statistical analyses with the Principal Component Analysis (PCA) and the Hierarchical Cluster Analysis (HCA) methods were performed on the complete chemical composition of the spontaneous volatile emissions of the six O. basilicum samples. The PCA was carried out on a 36 × 6 covariance data matrix (36 compounds × 6 samples = 216 data), selecting the two highest PCs obtained by the linear regressions, and the chosen PC1 and PC2 covered 61.30 and 25.90% of the variance, respectively, for a total explained variance of 87.20%. The HCA was performed using Ward's method, with squared Euclidean distances as a measure of similarity. Moreover, the statistically significant differences between content of the chemical classes in the analyzed samples were also evaluated with the analysis of variance (ANOVA) test. Averages were separated by Tukey's post hoc test and a p < 0.05 was used to assess the significance of differences between means.
The statistical analyses were performed using the JMP Pro 14.0.0 software package (SAS Institute, Cary, NC, USA).

Data Availability Statement:
The data that support the findings of this study are available on request from the first author Marina Laura. The DNA sequence presented in this study is openly available in GenBank, accession n. MT319764.